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هدف از این پژوهش ارزیابی و تفسیر اثر برهمکنش ژنوتیپ × محیط برای لاین (ژنوتیپ) های پیشرفته عدس ديم در 
هشت محبط آزمایشی بود. برای این منظور. تعداد ۱۵ لابن پیشرفته عدس ديم به‌همراه دو رقم شاهد (گچساران و سیهر) در 
WE‏ طرح بل وک‌های fol‏ تصادفی با سه تکرار در چهار ایستگاه تحقیقاتی گچساران, مغان» خرم‌آباد و ایلام در دو سال 
زراعی ۱۳۹۸-۹۹ و ۱۳۹۹-۱۴۰۰ بررسی و ارزیابی شدند. میانگین عملکرد دانه کل لاین‌ها ۸٩۱‏ کیلو گرم در هکتار بود. ژنوتیپ 
۰ با ۱۰۳۱ کیل وگرم در هکتار بیشترین و ژنوتیپ ۱۶ (رقم گچساران) با VIA‏ کیل و گرم در هکتار کمترین میانگین عملکرد را در 
مجموع Luxe‏ ها دارا بودند. با توجه به معنی‌دار بودن اثر برهمکنش ژنوتیپ × محیط. تجزیه پایداری عملکرد دانه برای 
ژنوتیپ ها pled!‏ شد. بررسی پایداری عملکرد دانه ژنوتیپ ها با روش رتبه (Rank)‏ شان داد که ژنوتیپ های ۱۲ و ۱۰ به 
تر تیب ژنوتیپ های با عملکرد دانه پایدار بودند. افزون بر این تجزیه پایداری با شاخص WAASB‏ به عنوان روش چند متغیره 
نیز انجام شد. بای پلات عملکرد دانه در برابر WAASB‏ نشان داد که اغلب ژنوتیپ ها دارای پایداری عملکرد دانه بیشتر از 
رقم های شاهد بودند. ژنوتیپ های ۱۲ ۰۱۳ ۱۱ ۱۰ وا به ترتیب دارای عملکرد دانه و پایداری عملکرد دانه بیشتر از رقم شاهد 
(سپهر) بودند. با در نظر گرفتن سهم برابر برای هر یک از دو جزء میانگین عملکرد و پایداری عملکرد دانه در محاسبه شاخص 
۷ زنوتیپ های ۱۰ ۱۳ ۸۱۱ ۱۲ ۱ به ترتیب در رتبه های برتر قرار گرفتند. 


واژه‌های کلیدی: عدس بای پلات» برهمکنش ژنو تیپ X‏ محبط. بهترین پیش‌بینی‌های ناار بب خطی. 323 4« AMMI‏ 
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از tee‏ روش‌های کاربردی در تفسیر 
برهمکنش ژنوتیپ × محیط در آزمایش‌های 
یکنواخت چند محیطی تجزیه (AMMI‏ 
Additive main effect and multiplicative‏ 
interaction)‏ است. این روش بواسطه امکان 
نمایش گرافیکی برهمکنش: نسبت به 
روش هایی همچون تجزیه واریانس ساده برتر 
است )2016 (Van Eeuwijk et al.,‏ تجزبه 
1 به طور گسترده ای در پژوهش‌های 
مربوط به ارزیابی برهمکنش ژنوتیپ × محیط و 
پایداری عملکرد دانه ژنوتیپ‌های گونه های 
زراعی گوناگون مورد استفاده قرار گرفته است 
(Rahayu, 2020: Hasani ef al., 2021;‏ 
Namdari et al., 2022a)‏ 

علاوه بر روش AMMI‏ روش بهترین 
پیش‌بینی‌های SENT‏ خطی (Best linear‏ 
unbiased perdiction = BLUP)‏ نیز برای 
برآورد میانگین عملکرد ژنوتیپ‌ها در 
آزمایش‌های چندمحیطی مورد استفاده می گیرد 
(Nardino et al., 2016; Olivoto et al.,‏ 
)2017 تجزیه AMMI‏ » بش عمده 
برهمکنش ژنوتیپ × محیط را در نخستین 
محور مولفه اصلی برهمکنش (1۳0۸1) توضیح 
می دهد. در حالیکه بیشترین خطای تصادفی 
به PCA‏ های پایانی مربوط می‌شود 
(Mofidian and Moghaddam, 2013)‏ 

برای برطرف کردن این نقطه ضعف روش 
AMMI‏ تلفق دوروش AMMI‏ و BLUP‏ 
پیشنهاد ote‏ است (Olivoto et al., 2019a)‏ به 


مقد مه 

مقدار و پراکنش بارندگی دوعامل عمده 
محدود کننده عملکرد als‏ عدس دیم در کشور 
می‌باشند و معرفی رقم های سا زگار با این 
محدودیت ها می تواند تاثیر به سزایی در بهبود 
عملکرد دانه و تولید عدس داشته باشد 
(Pezeshkpour and Afkar, 2019)‏ در زمينه 
uoces ills Glas pleads‏ 
معتدل در کشور از نظر تعداد رقم‌های معرفی 
شده محدودیت‌هایی وجود دارد. از جمله تا 
کنون رقم‌ه‌ای "گچساران» سپهر و 
بیله‌سوار "مناسب جهت کشت در دیم‌زارهای 
مناطق گرمسیری و نیمه گرمسیری معرفی 
شده‌اند. 

پایداری عملکرد دانه به مفهوم OLS‏ عملکرد 
دانه یک ژنوتیپ در سال های مختلف و 
محیط‌های گوناگون است و به این منظور 
در فرایند معرفی رقم (AA A‏ عملکرد دانه 
ژنوتیپ‌ها در محیط‌های گوناگون و در SB‏ 
آزمایش‌های چندمحبطی مورد ارزیابی قرار 
می گیرد» تا امکان سنجش برهمکنش ژنوتیپ × 
محیط فراهم شود. وجود برهمکنش بالای 
ژنوتیسپ X‏ محیط بویژه در مورد صفات 
پیچیده‌ای همچون عملکرد دانه» موجب 
کندشدن فرایند گزینش ژنوتیپ‌ها و معرفی 
رقم cla‏ جدید می‌شود )2020 (Sharifi,‏ 
بنابراین تفسیر درست از این برهمکنش در 
فرایند ارزیابی و شناسایی ژنوتیپ های برتره 
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بدهد. 

معرفی و استفاده از شاخص WAASB‏ برای 
ارزیابی پایداری عملکرد دانه ژوتیپ‌ها 
در آزمایش‌های چندمحیطی اخیرا مورد 
توجه به نژادی گران قرار گرفته است. از 
این رو تاکنون پژوهش‌های معدودی با 
استفاده از این شاخص برای بررسی 
برهمکنش ژنوتیپ X‏ محیط انجام شده است. 
در این‌باره اخیرا در پژوهش‌های به‌نژادی 
مرتبط بابرهمکنش ژنوتیپ × محیط در 
گندم دوروم )2021 (Karimizadeh et al.,‏ « 
عدس g(Namdari et al, 2022a)‏ 
برنج )2021 (Sharifi et al.,‏ از شاخص 
8 به منظور بررسی پای‌داری 
عملکرد دانه ژنوتیپ‌های آزمایشی استفاده 
شده است. 

هدف از ان پژوهش» شناسایی 
ژنوتیپ(های) برتر عدسدیم واجد عملکرد 
دانه و پایداری عملکرد دانه بالا برای 
معرفی رقم(های) جدید بود. بدین منظوره 
ارزیابی پایداری عملکرد RY‏ (ژنوتیپ) های 
پیشرفته عدس دیم با بهره گیری از روش تک 
متغیره رتبه (Rank)‏ و روش چند متغیره با 
محاسبه دو شاخص WAASB‏ و WAASBY‏ 


انجام شد. 


مواد و روش‌ها 
۱۳۹۸-۹ و ۱۳۹۹-۱۴۰۰ در چهار ایسستگاه 


wa 


همین منظور دو روش ذکرشده با معرفی 
شاخصی با نام "شاخص پایداری میانگین وزنی 
نمرات مطلق بهترین پیش‌بینی‌های نااریب خطی 
(Weighted average of absolute scores =‏ 
WAASB)‏ " تلفیق می‌شوند. برای محاسبه 
۵ به جای Sl‏ 4 تجزیه AMMI‏ روی 
داده‌ها در آزمایش چند محبطی انجام pags‏ 
روی بهترین پیش بینی های ااریب خطی 
(BLUPS)‏ انجام می شود. از نقطه نظسر 
کاربردی» AMMI BLUP‏ را می‌توان به 
عنوان دو رویکرد متمایز برای رسیدن به هدفی 
OLS‏ (متمایز نمودن الگوی برهم کنش 
ژئوتیپ ‏ محیط از خطای تصادفی) در نظر 
gta.‏ 

مزیت های شاخص WAASB‏ در مقاسه با 
شاخص‌های تجزبه oT AMMI‏ است که بر 
انحراف‌های Glee‏ )4 جای e^‏ انحراف‌ها در 
تجزیه (AMMI‏ استوار است. و حساسیت 
کمتری به داده‌های پرت دارد. افزون بر این 
در محاسبه WAASB‏ همه محورهای 
مولفه های اصلی مورد استفاده قرار می گیرند 
(Olivoto er al, 2019a)‏ پس از 
محاسبه شاخص WAASB‏ امکان محاسبه 
شاخص دیگری به نام WAASBY‏ نیز 
وجود دارد که در آن عملکرد دانه و پایداری 
عملکرد ژنوتیپ‌هاء همزمان در نظر گرفته 
می‌شوند. و با توجه به هدف برنامه به نژادی» 
به‌نژادگر می تواند وزن‌های متفاوتی به 
شاخص پایداری (WAASB)‏ و عملکرد(۷) 
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استفاده شد. در محاسبه شاخص میانگین وزنی 
نمرات مطلق بهترین پیش‌بینی‌های نااریب خطی 
(WAASB)‏ تجزبه AMMI‏ 45 جای al‏ بر 
روی داده های خام انجام resus‏ 
نحوی که دو روش AMMI‏ و L BLUP‏ 
یکدیگر تلفیق شدند. در این محاسبات از 
آماره‌های تجر 4 AMMI‏ 3 نیز آماره‌های oe‏ 
بر تجزیه AMMI‏ و ماتریس برهمکنش ژنوتیپ 
× محیط حاصل از بهترین پیش‌بینی‌های ااریب 
خطی (BLUP)‏ استفاده شد. 

نحوه محاسبه آماره‌های تجزیه AMMI‏ در سایر 
پژوهش‌ها تشریح (Karimizadeh et al, ilo‏ 
Olivoto et al., 20190: Pezeshkpour‏ ;2021 
et al, 2021 )‏ در اینجابه محاسه دو آماره 
WAASBY E WAASB‏ اشاره می‌شود: 


Lye PCA; XEP, | :۱ رابطه‎ 


WAASB;- - 
xı EP 


در این رابطه WAASBi‏ می‌انگین وزنی 
امتیاز های مطلق مربوط به ژئوتیپ شماره 1 است. 
ڄ1g ۴٣۵‏ امتیاز ژنوتیپ شماره 1 در محور مولفه 
اصلی برهمکنش شماره )می باشد. EPk‏ 
واربانس توجبه شده توسط IPCA‏ شماره ۸ است. 

به منظور فراهم آمدن امکان گزینش ژنوتیپ ها 
بر پایه هر دو شاخص عملکرد دانه و پایداری 
عملکرد» شاخص WAASBY‏ نیز محاسبه شد که 
فرمول OT‏ در رابطه ۲ ارائه شده است. 
رابطه UN‏ 


تحقیقات کشاورزی دیم گچساران؛ خرم‌آباد؛ 
مغان و ایلام در قالب طرح بل وک های کامل 
تصادفی با سه تکرار اجرا شد. مجموع تعداد 
محیط cls‏ آزمایشی هشت محیط شامل چهار 
مکان و دو سال بود. 

کد های محیط های آزمایشی به ترتیب: 
گچسار ان سال اول و دوم (G99 G98)‏ 
خرم‌آباد سال اول و دوم (K99 K98)‏ مغان 
سال اول و دوم «(M99 M98)‏ ایلام سال اول و 
دوم )198 199( بود. پانزده لاین (ژنوتیپ) 
پیشرفته پس از طی نمودن سه سال ارزیابی در 
قالب خزانه های بین المللی» آزمایش مقدماتی 
مقاسه عملکرد دانه» و آزمایش مقاسه عملکرد 
دانه پیشرفته» برای ادامه ارزیابی در قالب آزمایش 
یکنواخت سراسری مقایسه عملکرد دانه» انتخاب 
شدند. نام و شجره pV‏ (ژنوتیپ) cle‏ انتخابی و 
منشا آنها در جدول ۱ اراثه شده است. 

وضعیت دما و بارندگی در سال‌های اجرای 
آزمایش در جدول Y‏ ارائه شده است. در JUS‏ 
ژنوتیپ‌های برگزیده دو رقم "گچساران "و 
"سپهر" به عنوان رقم‌های شاهد (ارقام تجاری در 
مناطق یادشده) نیز در آزمایش مورد استفاده قرار 
گرفتند. کشت در همه ایستگاه‌های تحفیقاتی در 
تناوب گندم-عدس انجام شد و تاریخ کشت نیز 
بسته به منطقه» از آغاز آذر تا اواسط دی سته به 
زمان دریافت بارش های موثر بود. 

برای تجزیه آماری از بسته ویژه آزمایش‌های 
چندمحیطی با نام METAN‏ در بستر نرم افزار R‏ 


| 
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جدول۱- نام /شجره و منشاء aay‏ (ژنوتیپ) های عدس m>‏ 


120161. Name/pedigree and origin of rainfed lentil lines (genotypes) 
کد ژنوتیپ منشا نام /شجره کد ژنوتیپ‎ opel منشا نام‎ 
Genotype code Name/Pedigree Origin Genotype code Name/Pedigree Origin 

G1 2009S 96575-1 ICARDA G10 2009S 96501-5-X2007S61-ILL5883 x ILL6458 ICARDA 
G2 08S 40111-01 ICARDA G11 ILL 4605-2 ICARDA 
G3 ILL 5988 ICARDA G12 08S 41205-13-X20058140-ILL5883 x ILL7620 ICARDA 
G4 ILL 4605-1 ICARDA G13 2009S 96101-4 ICARDA 
G5 ILL 6183 ICARDA G14 X2012S-154-X2012S-154-ILL6002 x ILWL118 ICARDA 
G6 ILL 590 ICARDA G15 X2012S-102-X2012S-102-ALL7986 x ILWL 74) x (ILLA605 x ILL5677) ICARDA 
G7 2009S 96505-2-X2007S67-ILL6434 x ILL8199 ICARDA G16 Gachsaran IRAN 
G8 2009S 96518-1-X2007S107-ILL7940 x ILL5883 ICARDA G17 Sepehr IRAN 
G9 2009S 96501-1-X2007S61-ILL5883 x ILL6458 ICARDA 


جدول ۲- میانگین ماهانه دما و بارندگی ایستگاه های تحقیقات کشاورزی در سالهای زراعی ۱۳۹۸-۹۹ و ۱۳۹۹-۱۴۰۰ 


Table 2. Mean monthly temperature and rainfall of agricultural research stations in 2019-20 and 2020-21 cropping seasons 


ايلام مغان خرم آباد گچساران 
Gachsaran Korramabad Moghan Ilam‏ 
دما (درجه سانتی گراد) بارش Ge de)‏ دما (درجه سانتی گراد) بارش Ce de)‏ دما (درجه سانتی گراد) بارش Ge de)‏ دما (درجه سانتی گراد) بارش (میلی متر) 
Rainfall (mm) Temperature (°C) Rainfall (mm) Temperature (°C) Rainfall (mm) Temperature (°C) Rainfall (mm) Temperature (C) __‏ ماه Month‏ 
2020-21 2019-20 2020-21 2019-20 2020-21 2019-20 202021 2019-20 2020-231 2019-20 202021 2019-20 2020-231 2019-20 2020-21 2019-20 
22.7 24.8 0.0 15.0 18.5 19.1 29.2 10.5 22.9 22.9 94 12.4 26.3 27.8 24 0.0 مهر Oct.-Nov.‏ 
16.7 14.2 33.7 44.6 13.8 12.1 6.2 3.4 13.0 13.0 11.6 54.0 20.0 19.1 33.3 64.5 آبان Nov.-Dec.‏ 
100 94 127.8 134.4 6.2 7.0 25.9 274 77 7/1 154.0 112.0 14.5 13.6 401.5 111.0 آذر Dec.-Jan.‏ 
7.5 7.1 6.9 37.4 4.6 6.1 18.8 15.5 64 6.4 49 69.0 11.5 11 10.2 59.4 دی Jan.-Feb.‏ 
9.0 6.4 97.4 44.3 6.0 63 20.0 18.6 5.6 5.6 87.3 36.0 13.1 11.2 67.1 41 بهمن Feb.-Mar.‏ 
TI 15.7 15.2 189.0 27.8 11.4 114 19.5 54.4 10.2 53 283.1 32.1 11.8 11.8‏ 35.4 اسفند Mar.-Apr.‏ 
16.0 14.0 4.3 33.5 13.7 11.1 17.4 223 13.0 13.0 6.5 39.6 21.1 17.6 0 97.2 فروردین Apr.-May‏ 
23.5 20.2 0.1 11.3 19.7 17.1 15.4 43.0 18.8 18.8 2.0 12.0 26.7 24.1 0 27 اردبهشت May -Jun.‏ 
302.3 603.6 187.3 160.2 303.5 524 5915 435.5 مجموع Total‏ 
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ابتدا تجزیه واریانس ساده برای عملکرد دانه 
در هر یک از محیط‌های آزمایشی انجام OLS‏ 
و همگنی واریانس خطاهای آزمایشی با آزمون 
بارتلت تایید شد. با توجه به معنی دار بودن 
برهمکنش ژنوتیپ ‏ محیط در تجزیه واریانس 
م رکب (جدول ۳ تجزبه پایداری عملکرد 
دانه و ارزیابی برهمکنش ژنوتیپ X‏ محیط» 
انجام گردید. میانگین عملکرد دانه کل 
ژنوتیپ ها ۱ کیل و گرم در هکتار در کلیه 
محیط ها بود. ژنوتیپ ۱۰ با ۱۰۳۱ کیل و گرم در 
هکتار بیشترین و ژنوتیپ ۱۶ (رقم گچساران) با 
۸ کیل و گرم در هکتار کمترین میانگین 
عملکرد دانه را داشتند. 


در این رابطه :۷۷۸۵۸۹13۷۲ میانگین وزنی 
شاخص WAASB‏ و عملکرد (GY) ails‏ برای 
ژنوتیپب شماره WY d‏ وزن داده شده 4 عملکرد 
GYi cals‏ میانگین عملکرد ژنوتیپ شماره d‏ 
GYmax‏ بیشترین میانگین عملکرد اندازه گیری 
شده WS‏ وزن داده شده به پایداری» WAASBiI‏ 
مقدار شاخص WAASB‏ برای ژنوتیپ شماره d‏ 
WAASBnin‏ کمترین مقدار شاخص WAASB‏ 
اندازه گیری شده در ميان ژنوتیپ ها( که برابر 
است با بیشترین پایداری) می باشند. 

به منظور فراهم آوردن امکان مقایسه دقیق 
لاین ها (ژنوتیپ‌ها) در محیط های آزمایشی با 
رسم شکل‌های گوناگون در بستر نرم افزار R‏ 
برهمکنش ژنوتیپ × محیط بررسی و تفسیر شد. 


Table 3. Combined analysis of variance for seed yield of lentil genotypes environments 


ميانگین مربعات ‏ درجه آزادی 


S.O. V. منبع تغییرات‎ d.f. Mean Squares 
Environment (E) lo 7 10985517** 
Replication/E JUSSI محیط‎ 16 269470 
Genotype (G) ژنوتیپ‎ 15 146351** 
GxE زژنوتیپ × محبط‎ 112 102881** 
Error es 6 50802 


***: Significant at the 1% probability level. 


رنگ سبز روشن تغیبرات ناشی از برهمکنش 
ژنوتیپ × محیط و بخش‌های به رنگ سبز تیره 
تغییرات ناشی از ژنوتیپ را OUS‏ می‌دهند. به جز 
در مورد مولفه دوم» در مورد سار مولفه ها سهم 
تغییرات ناشی از برهمکنش ژنوتیپ × محیط 
فاصله بسیار زیادی با سهم اثر ژنوتیپ به تنهایی 
دارد (Fig. 1b)‏ 
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e‏ معنی دار در سطح احتمال ER‏ درصد. 


برهمکنش ژنوتیپ × محیط به هشت مولفه 
تقسیم شد که در شکل ۱ سهم هریک از مولفه‌ها 
در توضیح برهمکنش ژنوتیپ محیط نمایش داده 
شده است (Fig. 1a)‏ سهم مولفه‌ها از واربانس 
ژنوتیپ و برهمکنش ژنوتیپ × محیط در شکل ۱ 
نشان داده شده است (Fig. 1b)‏ اگر کل سطح 


مربع را ۱۰۰ درصد در نظر بگیریم» بخش‌های به 
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Fig. 1. (a) Contribution of each principle component in explaining GE interaction, and 


(b) contribution of principle components to variation of G+GE 


ایلام و خرم آباد مقدار بارش در دو سال زراعی 
تفاوت قابل توجهی داشت. که منجر به تفاوت 
قابل توجه در عملکرد دانه ژنوتیپ ها شد. 

در مجموع دو سال» در ایستگاه گچساران 
ژنوتیپ ۱۰ با میانگین عملکرد ۱۴۲۲ کیل وگرم 
در هکتار بیشترین عملکرد دانه را داشت» و با 
فاصله قابل توجه ای ژنوتیپ AM ۸۱۱ cle‏ ۴ و 
۲ تقریبا در سطح GUS‏ قرار داشتند (جدول 
۴ در ایستگاه خرم LT‏ ژنوتیپ های ۰۳ ۸۱۲ 
۳ و ابیشترین عملکرد دانه راداشتند. در 
ایستگاه مغان ژنوتیپ های ۸و ۱۲ و در ایستگاه 
ایلام بیشترین عملکرد دانه متعلق به ژنوتیپ های 
۰ ۱ بود (جدول ۴). بابد توجه داشت که 


برای شناسایی ژنوتیپ oh‏ برتر توجه به هر دو 
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نکته ای که در رابطه با عملکرد دانه 
ژنوتیپ های عدس در محیط ها مختلف باید 
مورد توجه قرار cola‏ این است که به سیب 
تفاوت های آب و هوایی قابل توجه بین 
cls Low‏ آزمایشی» عملکرد دانه ژنوتیپ‌ها در 
محیط‌های آزمایشی تغییرات قاببل توجهی 
داشت که نشان دهنده وجود اثر معنی دار 
برهمکنش ژنوتیپ X‏ محیط بود (جدول Y‏ 
برای مثال» در ایستگاه گچساران با وجود بارش 
بیشتر در سال زراعی دوم (جدول «Y‏ به Ho‏ 
پرا کنش امناسب و روی دادن قسمت عمده این 
بارش ها پیش از کاشت. در واقع دوره رشد 
کوتاه شد و گیاه با تنش خشکی قابل توجهی در 
انتهای فصل رشد مواجه گردید. در دو ایستگاه 
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ژنوتیپ های عدس در محیط های گوناگون را 
گزارش کرده اند :2006 (Bicer and Sarkar,‏ 
Kumar et al., 2005; Sabbaghnia et al.,‏ 
)2008. 


خصوصیت عملکرد ails‏ بالا و پایداری عملکرد 
لازم است. به همین دلیل» تجزیه پایداری 
عملکرد دانه امری ضروری است . سایر 


پژوهشگران نیز واکنش متفاوت عملکرد دانه 


جدول ۴- میانگین دو ساله عملکرد دانه (کیلو گرم در هکتار) ژنوتیپ های عدس در چهار مکان 


Table 4. Two-year mean seed yield (kg ha") of lentil genotypes in four locations 


v5 گچساران‎ the مغان‎ eu v5 ايلام مغان خرما باد گچساران‎ 
Genotype  Gachsaran  Koramabad  Moghan Iam Genotype  Gachsaran Korramabad Moghan Ilam 
Gl 1131 1302 467 817 G10 1432 1128 357 1206 
G2 1136 1121 447 736 Gli 1218 1200 418 1150 
G3 796 1378 538 936 612 1130 1308 566 1078 
G4 1174 1073 427 583 G13 1189 1310 439 863 
G5 900 1137 470 831 G14 1196 1149 428 831 
G6 903 1089 413 1052 G15 1083 1088 441 699 
G7 812 1250 451 866 G16 911 1032 413 559 
G8 795 1169 615 1008 G17 996 925 506 1090 
G9 831 1339 440 813 


LSD 5% = 115 kg ha" 


نداشت. پس از آن» ژنوتیپ های ESAGY‏ به 
عنوان ژنوتیپ chlo cls‏ عملکرد ails‏ پایدار 
Ab Sede) ch, alate‏ 

ناپایدارترین عملکرد دانه به ژنوتیپ های ۱۶ 
(رقم گچساران) و ۶ تعلق داشت که عمدتا 
در بخش پایینی رتبه بندی قرار داشتند 
(جدول ۵). موثر بودن آماره TOP‏ در 
شناسایی ژنوتیپ های پرمحصول وبا 
عملکرد پایدار توسط ple‏ پژوهشگران 
نیز گزارش شده است )2008 (Ebadi et al.,‏ 
در بسیاری از موارد. ژنوتیپ glia‏ با 
عملکرد پایدار از نظر آماره TOP‏ از نظر 
عملکرد نیز در رتبه (SVG cls‏ قرار داشتند 
(Abdullahi et al., 2007)‏ 
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تجزبه پابداری عملکرد دانه با استفاده از روش 
ناپار امتری 

آماره ناپارامتری فاکس )1990 (Fox et al.,‏ 
ss)‏ ها براساس رتبه دریافتی در سه رده 
(Top) Yu‏ پایین (Mid) «tle (Low)‏ در 
محیط cle‏ آزمایشی رتبه GAL‏ می کند و بر 
مبنای درصد قرا رگیری ژنوتیپ ها در سه بخش» 
ژنوتیپ با پایدارترین عملکرد شناسایی می شود. 
بالاتر بودن مقدار آماره TOP‏ یک ژنوتیپ به 
معنای پایدارتر بودن عملکرد دانه آن ژنوتیپ 
است. بر مبنای آماره TOP‏ ژنوتیپ Y‏ دارای 
پایدارترین عملکرد دانه بود» چون در VO‏ درصد 
موارد در رده بالایی از نظر رتبه بندی قرار 


داشت» و هرگز در بخش پایینی رتبه بندی قرار 
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جدول ۵- رتبه بندی عملکرد ژنوتیپ های عدس با استفاده از روش فا کس و همکاران (۱۹۹۰) 


Table 5. Ranking of seed yield of lentil genotypes using Fox et. al. (1990) method 


بالا wh se‏ ژنوتیپ بالا owh le‏ ژنوتیپ 

Genotype Low Mid Top Genotype Low Mid Top 
G1 50.0 25.0 25.0 G10 25.0 25.0 50.0 
G2 50.0 25.0 25.0 G11 25.0 37.5 37.5 
G3 25.0 50.0 25.0 G12 0.0 25.0 75.0 
G4 50.0 125 37.5 G13 12.5 50.0 37.5 
G5 50.0 50.0 0.0 G14 375 25.0 37.5 
G6 62.5 12.5 25.0 G15 37.5 62.5 0.0 
G7 37.5 50.0 12.5 G16 62.5 37.5 0.0 
G8 12.5 37.5 500 G17 12.5 62.5 25.0 
G9 50.0 125 37.5 


وزن هایی به این La Sl‏ داده می‌شود. در OLALS‏ 
زراعی مانند ذرت )2016 (Baretta et al.,‏ و 
eslex_uI« (Barbosa ef al., 2014) ye‏ از 
روش‌های مبتنی بر BLUP‏ در ارزیابی کارایی 
ژنوتیپ ها و پیش بینی دقیق اجزای واریانس 
ژنوتیپی مژثر بوده اند. 

به منظور ارزیابی پایداری عملکرد دانه 
ژنوتیپ cle‏ عدس از چند شکل بهره گرفته شد 
که یکی از آنها انحراف عملکرد دانه اسمی هر 
یک از ژنوتیپ‌ها با عملکرد پیش‌بینی شده در 
روش BLUP‏ بود (شکل Y‏ در این «JS—‏ 
خطهایی که با رنگهای گوناگون رسم شده اند 
مربوط به هریکک از ژنوتیپ ها می باشند» و 
هرچه شیب خط بیشتر باشد OUS‏ دهنده انحراف 
بیشتر از عملکرد دانه پیش بینی شده و پایداری 
عملکرد کمتر می باشد. این بای‌پلات نمایش 
دیگری از بای‌پلات "کدام ژنوتیپ- کجا 
برنده " در روش چند متغیسره دیگر تجزیسه 
پایداری GGE)‏ بای پلات) می باشد. 


تجزبه پایداری بااستفاده از آماره SLA‏ چند 
WAASB o jio‏ و WAASBY‏ 

عملکرد دانه برآورد شده برای هر یک از 
زنوتیپ‌های عدس باروش BLUP‏ در شکل۲ 
نمایش داده شده است. برآورد دقیق میانگین 
عملکرد دانه ژنوتیپ ها بویژه در مدل های مختاط 
از جمله مزیت cle‏ عمده روش  BLUP‏ باشد 
(Olivoto et al., 2019a)‏ افزون بر این در 
صورتی که اثر مختلط خطی وجود داشته باشد» 
BLUP‏ امکان پیش بینی ang‏ اثرهای تصادفی را 
امکان پذیر می سازد )2005 (Smith er al.,‏ 
روش BLUP‏ کارایی خود را در بهبود دقت پیش 
بینی اثرهای تصادفی ابت e» S.‏ است (Piepho et‏ 
al. 2008)‏ روش BLUP‏ این ویژگی را دارد که 
میانگین ها را با دقت YL‏ به‌ویژه در مدل‌های 
مختلط بر آورد می کند ;2016 (Nardino et al.,‏ 
Olivoto et al., 20199: Olivoto et al.,‏ 
(20192 در روش BLUP‏ در وهله نخست 


اثرهای مدل ANOVA‏ بر آورد می‌شوند» و سپس 
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شکل ۲- بهترین پیش بینی های نااریب خطی برای W‏ ژنوتیپ عدس. دایره های آبی و قرمز به ترتیب 
نشان دهنده ژنوتیپ هایی با ميانگین عملکرد دانه بالاتر یا پایین تر از میانگین های عملکرد دانه پیش بینی 
شده با BLUP‏ می باشند. 

Fig. 2. Best linear unbiased prediction (BLUP) for 17 lentil genotypes. Blue and red 


circles indicate above and below mean seed yield than BLUP predictions, respectively. 
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نخسن مولفه اصلی محیطی 
Environment 1‏ 


(IPCAI) محبط‎ X اصلی محبطی برهمکنش ژنوتیپ‎ ail ge شکل ۳- عملکرد دانه اسمی در برابرامتباز های نخستین‎ 
Fig. 3. Nominal seed yield vs. first principal component environmental axis of genotype 


x environment interaction (IPCA1 ) 
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سا زگاری خصوصی ژنوتیپ ها استفاده شده 
(Gurmu et al., 2012; Yan et al, Cul‏ 
Farshadfar, 2013)‏ ;2000 

محاسبه شاخص WAASBY‏ با در نظر 
گرفتن سهم های گوناگون برای هر یک از 
دو جزء عملکرد دانه و پایداری عملکرد 
ژنوتیپ های عدس انجام شد AF JS)‏ 
همانگونه که درشکل ۴ مشاهده می شود 
رتبه بندی ژنوتیپ ها با شاخص WAASBY‏ 
بسته به اينکه چه نسبتی از شاخص WAASB‏ و 
عملکرد دانه در نظر گرفته شود تفاوت زبادی 
پیدا می کند. هر چه از سمت چپ به سمت 
راست > OS‏ کنیم سهم پایداری عملکرد دانه 
در واحدهای پنج درصدی کاهش می bh‏ و به 
سهم عملکرد دانه افزوده می شود. 

با در نظرگرفتن سهم ۱۰۰ درصدی برای 
EEE‏ و di‏ رت کم 
عملکرد دانه همانگونه که در انتهای سمت 
چپ شکل ۴ آشکار است. ژنوتیپ‌های ۱۳ P‏ 
۵ و به ترتیب پایداری عملکرد دانه بالاتری 
داشتند. در SUE‏ با در نظر گرفتن عملکرد دانه 
s‏ تنھایی (انتهای سمت راست شکل 4۳ 
ژنوتیپ های ۱۰۰۱۲ و ۱۱بالاترین رتبه ها 
را به خود اختصاص دادند. با در نظر گرفتن 
سهم برابر (۵۰/۵۰) برای پایداری عملکرد و 
عملکرد دانه» در میانه شکل دیده می‌شود که 
ژنوتیپ های ۰۱۰ V OY ۰۱۱ OY‏ در رتبه های 
برتر قرار گرفتند. یکی از مزیت های 
شاخص WAASBY‏ این است که بسته به هدف 
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شماره های روی محور افقی از یک تا هشت 
نشان دهنده شماره هر یک از محیط های 
آزمایشی می‌باشند. و بر این اساس امکان 
ارزیابی عملکرد دانه هر یک از ژنوتیپ ها در 
هر محیط وجود دارد (شکل ۳. این ویژگی 
بویژه در مورد شناسایی وتعیین سا زگاری 
خصوصی زنوتیپ‌هامورد استفاده قرار 
می گیرد. با در نظر گرفتن موارد اشاره 
شده می‌توان چنین استنباط کرد که در میان 
زنوتیپ‌های عدس. ژنوتیپ های A‏ ۰۱۷ ۱۶ ۴ 
ود deu‏ رین نت E‏ بو ده و سای انم 
بیشترین تاثیر را در برهمکنش ژنوتیپ محیط 
داشتند. ژنوتیپ های ۱۶ و ۱۷ به ترتیب مربوط 
به رقم cle‏ گچساران و سپهر می باشند. 

ژنوتیپ هایی مطلوب شناخته می‌شوند که 
افزون بر داشتن عملکرد بالا (دستکم بالاتر از رقم 
های شاهد) پایداری عملکرد دانه بالایی نیز در 
cl Low‏ گوناگون داشته باشند و شیب خط 
ترسیم شده برای آنها کم باشد. در نقطه مقابل دو 
ژنوتیپ ۱۰ و ۱۱ کمترین شیب خط و بیشترین 
پایداری عملکرد ale‏ در محیط های آزمایشی 
داشتند. ply‏ ژنوتیپ ها نیز در محدوده متوسط 
قرار داشتند. همانگونه که در محور افنقی شکل 
اشاره شده است. ترسیم نمودار بر پایه مولفه 
ال نخس ی Obl‏ و سایر Ol js baci x‏ 
منظور نشده اند بنابراین تفسیر این شکل را باید با 
دقت پیشتری انجام داد. )2021 (Sharifi et al.,‏ 
از نمودار "کدام ژنوتیپ -کجا برنده " در 


پژوهش cle‏ گوناگون به منظور ارزیابی 
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ناپایدارترین عملکرد ails‏ را داشتند. 

در بخش دوم ژئوتیپ هاو محیط های با 
عملکرد دانه بالاتر از میانگین و پایداری 
عملکردکمتر قرار می گیرند. در این بخش 
ژنوتیپ LA‏ فاصله ای بسیار اند ک از میانگین 
عملکرد دانه و شاخص پایداری کمتر و ژنوتیپ 
۰ با میانگین عملکرد دانه بالا و شاخص 
WAASB‏ کمی بیشتر از میانگین (پایداری 
عملکرد کمتر از میانگین) قرار گرفتند. در مورد 
ژنوتیپ ۱۰ بايد اشاره شود که این ژنوتیپ از نظر 
عملکرد دانه در کنار ژنوتیپ MY‏ در رتبه نخست 
از نظر عملکرد دانه قرار گرفت و با توجه به 
فاصله بسیار کم از میانگین شاخص WAASB‏ 
برتری چشمگیری در مقایسه با هر دو رقم شاهد 
(ژنوتیپ cla‏ ۱۶ و ۱۷) داشت (شکل ۵. 
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برنامه به‌نژادی» می توان در محاسبه OT‏ به 
هریک از دو جزء عملکرد و پایداری عملکرد 
وزن‌های متفاوت داد تا سهم هر جزء پررنگ‌تر 
یا کم رنگ‌تر شود (شکل AF‏ 

بای پالات شاخص پابداری WAASB‏ در 
برابر عملکرد دانه در شکل ۵ نمایش داده شده 
است. همانگونه که در شکل ۵ مشاهده می شود 
این شکل به چهار بخش تقسیم شده است. در 
بخش ۱ ژنوتیپ ها و محیط هایی با عملکرد 
NER‏ وار فت ای 
8 بالا قرار دارند. مقدار بالاتر 
olis 8‏ دهنده پایداری عملکرد کمتر 
است. ژنوتیپ های ۰۱۷ ۱۶» ۴و ۱۵ در این 
بخش قرار گرفتند (شکل ۵). دو رقم شاهد 
گچساران و سپهر (ژنوتیپ های ۱۶ و OV‏ تقریبا 
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WAASB/SY نبت‎ 


WAASB/SY 


شکل ۴- رتبه بندی ۱۷ ژنوتیپ عدس با در نظ ر گرفتن نسبت های متفاوت شاخص پایداری عملکرد 


(SY) ails و عملکرد‎ (WAASB) 


Fig. 4. Ranking of lentil genotypes based on different stability index (WAASB) and 


seed yield (SY) 
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شکل ۵- بای پلات عملکرد ala‏ × شاخص پایداری WAASB‏ برای ۱۷ ژنوتیپ عدس در هشت محیط 


(چهار مکان و دو سال) 


Fig. 5. Seed yield (SY) x stability index (WAASB) biplot for 17 lentil genotypes in 


eight environments (Four locations and two years) 


بخش قرار گرفته اند. ژنوتیپ های ۱۴ و 3€ 
نظر عملکرد بسیار نزدیک به میانگین هستند و از 
این نظر با رقم شاهد سپهر تقریبا در رتبه برابری 
قرار داشتند. بر این اساس و با توجه به اهمیت 
عملکرد دانه» این ژنوتیپ ها از مطلوبیت کمتری 
od dd s uia‏ ها M AT AY‏ 5 383 
این بخش از بای پلات مطلوب ترند (شکل (à‏ 

بای پلات SY x WAASB‏ تا حدی مشابه 
بای پلات شناخته شده AMMIT‏ است؛ اما در 
مقایسه با OT‏ واجد مزیت هایی می باشد. چون 
در ترسیم OT‏ به جای اینکه تنها از مولفه اصلی 


Ww 


همچنین می توان ژنوتیپ ۱۰ را در کنار چند 
ژنوتیپ Ko‏ به عنوان یکی از ژنوتیپ های برتر 
و مناسب برای نامگذاری و آزاد سازی به Ol ge‏ 
رقم جدید در نظر گرفت. 

در بخش سوم ژنوتیپ ها و محیط هایی با 
عملکرد پایین و پایداری عملکرد بالا قرار 
دارند» که شامل ژنوتیپ A cle‏ ۲۶۷و b‏ 
می شوند. بخش چهارم مطلوب ترین ژنوتیپ ها 
را شامل می شود که همزمان عملکرد دانه و 
پایداری عملکرد بالاتر از ميانگین دارند (شکل 


۵ ژنوتیپ های ۰۱۲ ۱۳ ۱۱ ۱و ۳و ۱۴ در این 
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هر دو روش می شود و بافته cle‏ حاصل از 


آزمایش‌های چند محیطی دقیق تر تفسیر 


می شوند )2021 (Ajay and Singh,‏ 
در رابطه با اهمیت در نظر گرفتن همه 
مولفه های برهمکتش گزارش شده است که 
رتبه بندی ژنوتیپ های پولاف هنگامی که تنها 
دو مولفه در نظر گرفته شدند با هنگامی که همه 
مولفه ها در ارزیابی مدنظر قرار گرفتند به 
شدت تغییر کرد به نحوی که ژنوتیپی که با در 

نظ رگرفتن دو مولفه از نظر پایداری عملکرد د 
رتبه نخست قرار داشت در ادامه به رتبه ششم 
سقوط کرد (Olivoto et al., 2019a)‏ سایر 
پژوهشگران نیز ميزان توجیه برهمکنش ژنوتیپ 
x‏ محیط با استفاده از مولفه نخست. پایین بوده 
است. از OT‏ جمله می توان به نیشکر با ۳۲/۹ 
درصد )2011 &)3«(Ramburan ef al.,‏ با 
۴ درصد )2009 (Balestre et al.,‏ « گندم 
با ۳۲/۱ درصد )2017 (Tigabu et al.,‏ و ۳۳/۱ 
درصد )2019 (Veenstra et al.,‏ اشاره کرد. 
در اینگونه موارد بدیهی است که تفسیر 
برهمکنش با در نظر گرفتن یک مولفه بایستی با 


احتباط v‏ بیشتری انجام گیرد و در نظر گرفتن 
تعداد ia Avia M‏ ی 
همانگونه که ay‏ پیشتر هم توضیح داده شد. با نتایج 


دقیق تری همراه است. 

در پژوهش حاضر مولفه اصلی نخست 
توجیه کننده بیش از ۴۰ درصد از واریانس 
برهمکنش ژنوتیپ X‏ محیط بود که میزان قابل 
قبولی است. بدیهی است که با در نظر گرفتن 


۱۳۲ 


نخست IPCAT)‏ استفاده شود از همه 
IPCA‏ ها استفاده می‌شود. بنابراین به ساير 
مولفه ها نیز امکان بروز در بای پلات را 
می دهد که به نوبه خود موجب تفسیر دقیق تر 
نتایج می شود. در مقایسه با بای پلات elas"‏ 
برنده -کجا" که در شکل ۴ نمایش داده شده 
است نیز این مزیت صادق میباشد» چون 
در OT‏ شکل نیز همانگونه که اشاره شد 
تنها مولفه اصلی نخست مورد استفاده قرار 
گرفته است. از آنجاییکه در محاسبه 
شاخص. ۷۸۸58 از همه مولفه cls‏ اصلی 
بهره برده می‌شود؛ برهم کنش ژنوتیپ X‏ محیط» 
دقیق تر uli‏ داده می شود (Olivoto et al.,‏ 
2019b)‏ پیش از این نیز از بای SY x c»,‏ 
۶ در ارزیابی پایداری عملکرد دانه 
ژنوتیپ cla‏ عدس به طور موثری استفاده شده 
است (Namdari et al., 2022a)‏ 

در مواردی که مولفه اصلی نخست. درصد 
نسبتاً کمی از برهمکنش ژنوتیپ X‏ محیط را 
توجیه می کند. تفسیر نتایج بای پلات AMMII‏ 
ممکن است با دقت کمتری همراه باشد» زیرا در 
این حالت بخش cla‏ پیچیده تر برهمکنش در 
سایر مولفه ها قرار می گیرد و OLS‏ ارزیابی 
آنها در | 
این نقطه ضعف. بای پلات SY x WAASB‏ با 


ین بای پلات فراهم نیست ,در پاسخ به 


در نظر گرفتن هر دو بخش ساده و پیچید 
از برهمکنش ژئوتیپ X‏ محیط امکان دستیابی 
به تفسیر بهتر را فراهم می کند. ادغام دو روش 
AMMI s BLUP‏ موجب افزایش کارایی 
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(Sharifi et al., 2021) مى شود‎ 


سپاسکزاری 
نگارند گان بدینوسیله از مسدیریت 
استگاه های تحشقات کشاورژی گساران؛ 


مغان» خرم آباد و PIL!‏ برای فراهم آوردن 


تعداد بیشتری از مولفه هاء نتایج دقیق تر وتفسیر 
آنها نیز با قابلیت اطمینان بیشتری همراه است. 
در این رابطه گزارش شده است که هرچه 
در تجزیه AMMI‏ از مولفه‌های اصلی بیشتری 
استفاده شود نتایج به دست آمده با 


به نوبه خود از برتری ارزیابی با شاخص 
WAASBY‏ در مقایسه با AMMI‏ ناشی 


امکانات لازم برای اجرای این پسژوهش 
مییاسگزازی مین کد 
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